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Abstract 

The Kohn variational principle and the (correlated) Hyperspherical Har- 
monics technique are applied to study the n — 3 H and p — 3 He scattering 
at zero energy. Predictions for the singlet and triplet scattering lengths are 
obtained for non-relativistic nuclear Hamiltonians including two- and three- 
body potentials. The calculated n — 3 H total cross section agrees well with the 
measured value, while some small discrepancy is found for the coherent scat- 
tering length. For the p — 3 He channel, the calculated scattering lengths are 
in reasonable agreement with the values extrapolated from the measurements 

made above 1 MeV. 
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In the last few years the scattering of nucleons by deuterons has been the subject of a 
large number of investigations. This scattering problem is in fact a very useful tool for testing 
the accuracy of our present knowledge of the nucleon-nucleon (NN) and three nucleon (3N) 
interactions. Noticeable progress has been achieved, but a number of relevant disagreements 
between theoretical predictions and experimental results remains to be solved 

It is therefore of interest to extend the above mentioned analysis to four nucleon scatter- 
ing processes. In this case, an important goal for both theoretical and experimental analysis 
is to reach a precision comparable to that achieved in the N — d case. This is particularly 
challenging from the theoretical point of view, since the study of A = 4 systems is notice- 
ably more complicated than the A = 3 one. Recently, accurate calculations of the alpha 
particle binding energy S 4 have been achieved 0-0. It has been shown that, with NN+3N 
potential models fitting the 3 H binding energy, no four-nucleon potential seems necessary 
to reproduce the experimental value of ||. Therefore, it is expected that NN and 3N 
interactions should be sufficient to describe the four nucleon scattering processes too. Thus, 
discrepancies between theory and experiment would be useful to gain further information 
on the nuclear interaction. For example, the polarization observables in the reaction p — 3 H 
are believed to be very sensitive to the spin-orbit interactions JJ. Moreover, four nucleon 
reactions play an important role also in Astrophysics and other subfields of physics. 

In this letter, the problem of n — 3 H and p — 3 He zero energy scattering is studied. The 
aim is to obtain accurate estimates of the corresponding scattering observables by using NN 
and 3N realistic interactions. The relevant quantities in n — 3 H zero-energy scattering are 
the singlet a s and triplet a t scattering lengths. They can be obtained from the experimental 
values of the total cross section <tt and the coherent scattering length a c , 
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a T = n(\a s \ 2 + 3\a t \ 2 ) , a c =-a s + -a t . (1) 

The n — 3 H cross section has been accurately measured over a wide energy range and 
the extrapolation to zero energy does not present any problems. The value obtained is 
or = 1-70 ± 0.03 b 0. The coherent scattering length has been measured by neutron- 
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interferometry techniques. The most recent values reported in the literature have been 
obtained by the same group; they are a c = 3.82 ± 0.07 fm || and a c = 3.59 ± 0.02 fm ||, 
the latter value being obtained with a more advanced experimental arrangement. Recently, 
the estimation of a c = 3.607 ± 0.017 fm has been obtained from p — 3 He data by using an 
approximate Coulomb-corrected R-matrix theory JIT]]. 



The corresponding quantities for p — 3 He scattering are more difficult to evaluate. Ap- 
proximate values have been determined from effective range extrapolations to zero energy 
of data taken above 1 MeV, and therefore suffer large uncertainties |ll|,[l2l| . 

From the theoretical point of view, the problem of the scattering of four nucleons has been 
considered for a long time (see ref. jl3) and references cited therein). The most widely used 



techniques are based on the Faddeev-Yakubovsky (FY) equations |14|-[T6|| and the Kohn- 



Hulthen variational principles [17|. In the latter case, the Resonating Group Method (RGM) 



has been used to parametrize the wave function (WF) fL8|,19|, but also the expansion of the 



WF on a Hyperspherical Harmonic (HH) basis has been investigated P0| . Calculations 
using the FY and HH techniques, which allow for the full description of the four body 
dynamics, were performed by using simple central or separable potentials. Only recently, 
the FY equations have been solved by adopting realistic NN potentials [f2T] . 

In the present paper, the wave functions of the scattering states are expanded in terms of 
the correlated Hyperspherical Harmonic (CHH) basis [[22|] and the Kohn-Hulthen variational 
principles are applied. Such a technique have been successfully used in the study of the N — d 
scattering below and above the deuteron breakup threshold. The present calculations follow 
exactly the same line followed in the N — d case described in ref. | 23| . Let us consider the 



p — 3 He scattering; the case of n — 3 H scattering can be easily obtained in the limit e 2 — > 0, 



where e is the unit charge (see also ref. pi|). The WF with total angular momentum J, 



parity n and total isospin T, T z can be written as 

where the index 7 denotes hereafter the set of quantum numbers J, n, T, T z . The first term 



of eq. (0) must be sufficiently flexible to guarantee a detailed description of the "core" 
of the system, when all the particles are close to each other and the mutual interaction is 
large; ^ goes to zero when the p — 3 He distance r p increases. It has been expanded in 
terms of CHH basis functions, following the procedure discussed in detail in ref . |j| . 

The second term $2s describes the asymptotic configuration of the system, for large r p 
values, where the nuclear p — 3 He interaction is negligible. The quantum number L is the 
relative orbital angular momentum; S is the spin obtained by coupling the spin 1/2 of 3 He 
to the spin of the fourth nucleon. The angular momenta L and S are coupled to give the 
total angular momentum J. In the present study the total isospin is T = 1. The function 
&1 S must be the solution of the two-particle Schroedinger equation appropriate for large r p 
values. It is convenient to introduce the following surface functions 

n[% = E {Y L {*i) [*jWds} JJt fe^]TT^t A) (r«) , (3) 
i=i 

where the product $j ke x £^ is the WF ^/3 He of the 3 He bound state (in the case of n — 3 H 
scattering, it is the WF \&3 H of the 3 H bound state). They are normalized to unity and are 
antisymmetrical for the exchange of any pair of particles j, k and £. Both \l/3 He and \I/3 H 
have been determined as discussed in ref. [^3J by using the CHH expansion for a three-body 
system. Within this scheme the WF and the binding energy B 3 are determined with high 
accuracy. For example, the B 3 evaluated for the different potential models considered in this 
paper agree within a few keV to the corresponding results obtained by solving the Faddeev 



equations [25.26 



In eq. (||), the spin (isospin) function of the unbound nucleon i is denoted by Xi (&)■ 
Moreover, is the distance between nucleon i and the center of mass of 3 He. The functions 
"Rx\ r i) °f ec i- © can ^ e t & ken to be the regular (A = R) and irregular (A = I) radial 
solutions of the two-body Schroedinger equation without nuclear interaction. They are 
analogous to those used in iV — d scattering [f23j . 

With the above definitions, the asymptotic WF is written as 
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*Zs = ti L % + £ , (4) 

L'S' 

where the matrix element 7 -Rff/ gives the relative weight between the regular and the ir- 
regular L'S' components. These elements are the reactance matrix (R- matrix) elements, 



except for some numerical factors ||23|| . The eigenvalues of the i?-matrix are tan ^5, where 
6ls are the eigenphase shifts of the 2S+1 Lj nTTz wave. 

The convergence of the expansion of the internal part is conveniently studied by 
grouping the functions of the basis in "channels" (a given channel contains CHH states with 
the same angular-spin-isospin quantum numbers). It is very useful to consider first the 
channels with orbital angular momentum values as low as possible. One channel at a time is 
included in the expansion of the number of the CHH functions belonging to that channel 
is increased until convergence is reached. If the contribution of that particular channel 
is found to be sizeable, the corresponding CHH functions are retained in the expansion; 
otherwise, they are rejected. Then, others channels are added and the convergence studied 
in terms of the total number of channels N c . This procedure results to be effective since i) 
the value of N c can be kept rather low, and %%) a small number of CHH functions is sufficient, 
except for few channels. In particular, for the states (S-wave, T = 1) considered here, the 
number of channels included finally in the wave functions is rather small (iV c « 6-^8). 
This is due mainly to the Pauli principle which prevents the overlap of the four nucleons. 
As a consequence, the internal part is rather small and does not require a large number of 
channels. 

The quantities to be determined in the WF (0) are the hyperradial functions entering 
the HH expansion of the internal part ^/J, and the matrix elements 7 -Rjj/'. For these, the 
Kohn or the Hulthen variational principles have been used. The Kohn variational principle 
establishes that the following functionals 

VR s d\ = - ^i&lvin - E \*ls) . ( 5 ) 

where 7 -Rff/ are the trial parameters entering eq. (§), must be stationary with respect to 
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variations of all the trial parameters of the WF. In eq. (|| E is the total (cm.) energy and 
M the nucleon mass. 

The form of the equations then derived and the procedure to solve them is completely 



analogous to those of ref. [23] and is not repeated here. With the Hulthen variational 



principle the asymptotic function is written in the form 

^LS — + U LL'* L L'S'~f > \°) 

L'S' 

where U = R~ l . The Kohn and Hulthen variational principles lead to essentially different 
equations. Therefore, if the solutions in the two cases turn out to be close to each other, we 
are quite confident that they are close to the true solution. 

The results for the singlet and triplet scattering lengths for n — 3 H scattering are given in 
table pj as a function of the number of channels included in the WF. The potential adopted in 
this case is the AV14 interaction [p^j , so that a direct comparison with the results obtained 
in ref. [FITJ by solving the FY equations can be made. From an inspection of the table, 
the rapid convergence with N c is evident; this fact reflects that i) the scattering lengths are 
mainly determined by the asymptotic part and ii) the CHH expansion basis is very effective. 
Moreover, there is a strict agreement between the converged values of the scattering lengths 
obtained by means of the Kohn and the Hulthen variational principles. Both estimates 



compare very well with the FY results of ref. ||21|| , which is a strong signal of the good 
accuracy of both calculations. 

The calculated singlet and triplet n — 3 H scattering lengths corresponding to different 
potentials models are plotted versus the corresponding 3 H binding energy in fig. |l|. The 



most recent experimental values P,|T0! of a s and a t have also been reported. The models 



including only NN forces are the AV14 [[27], AV8 |28|] and AV18 |29|] potentials. Including 
3N forces we have : the AV14+Urbana model VIII (AV14UR) 0, AV18+Urbana model 
IX (AV18UR) 0, AV14+Brazil with A = 5.6™^ (AV14BR1) and AV14+Brazil with A = 
5.8?/^ (AV14BR2) 0. In the AV14UR and AV18UR models, one the parameters of the 
3N potentials was chosen so that to reproduce the experimental 3 H binding energy value 



B 3 = 8.48 MeV. The AV14BR1 and AV14BR2 models have been chosen so as to give slightly 
larger B 3 values. It should be noted that all the results for the singlet (triplet) scattering 
length fall essentially on a straight line. However, the experimental values extracted from 
the data do not lie on the theoretical curves. This disagreement is related to a rather small 
discrepancy between the calculated and measured coherent scattering length, as will be 
shown below. 

The calculated total cross section and coherent scattering length for the AV14UR and 
AV18UR models are compared with the experimental values |7HlO| in table |T[ These two po- 
tential models are chosen since they well reproduce the experimental B% value, and meaning- 
ful comparisons with the scattering data extracted from experiments can be then performed. 
From inspection of table ||, it can be concluded that there is a satisfactory agreement be- 
tween the calculated and the measured value of <jt- The calculated coherent scattering 
lengths, differ, however, by about 3% from the experimental values. This small discrepancy 
gives rise to the large differences in the scattering lengths, when these are determined from 
the relations given in eq. (|I]). In fact, in the a s , a t plane, the ellipse corresponding to the 
experimental values of the total cross section a? = 1.7 b and the straight line corresponding 
to the coherent scattering length a c = 3.7 fm are almost tangent. Therefore, a slight change 
in the a c value produces a large variation of a s and at. This is also the reason for the large 
uncertainty in the values of a s reported in figure [3]. 

The 3 He binding energy _B 3 ( 3 He) and the p — 3 He scattering lengths as determined with 



the AV18 and AV18UR models are presented in table [TTT| , together with the available ex- 
perimental data |IT| , [r2| . It should be remarked that, in contrast with the AV18UR model, 
the AV18 potential does not reproduce correctly the experimental value of I?3( 3 He). More 
in general, it has been verified that the scattering length values show a scaling property 
analogous to that found in the n — 3 H case. In table [TTI|, the available experimental values 
have also been reported. However, it should be observed that i) such experimental values 
have been extrapolated to zero energy from measured data taken above 1 MeV; ii) the 
quoted "error bars" include only statistical and not systematical uncertainties [32|. The 



p — 3 He experimental scattering lengths therefore suffer large uncertainties, even somewhat 
bigger than those reported in the table. By inspection of the table it can be concluded that 
the agreement between the AV18UR predictions and the experimental values is reasonably 
satisfactory and that it would be very useful to have a more accurate experimental deter- 
mination of a s and a t . Finally, it should be noted that the p — 3 He scattering lengths are 
larger than the corresponding values found in the n — 3 H case. This result is quite similar 
to that found in s-wave N — d scattering in the quartet spin state. 

In conclusion, accurate predictions of the n — 3 H and p — 3 He zero energy scattering 
lengths with realistic hamiltonians including NN and 3N potentials have been produced. The 
Kohn-Hulthen variational principle and the correlated Hyperspherical Harmonics technique 
were used to solve the four-body problem and to calculate the quantities of interest. The 
singlet and triplet scattering lengths for n — 3 H scattering were found to lie on straight 
lines when plotted against the 3 H binding energy for a variety of potential models. Our 
total cross section agrees well with the measured value, while some discrepancy is found 
in the comparison of the coherent scattering length values quoted in the literature. This is 
somewhat surprising, since the corresponding quantity in N — d scattering is well reproduced 
by the theory |2E], and the same was expected for the four-nucleon case. 

Although low-energy p — 3 He and n — 3 H experiments are difficult, we hope that the 
present work might inspire further efforts in this area. 

The authors wish to thank J. Carbonell and L.D. Knutson for valuable discussions and 
for providing their results prior to publication, and I. Bombaci and L. Lovitch for the critical 
reading of the manuscript. 
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FIGURES 

FIG. 1. Singlet (full symbols) and triplet (open symbols) scattering lengths plotted against the 
3 H binding energy. Circles labelled by a, b, c, d, e, f correspond to the AV18, AV14, AV8, AV18UR, 
AV14BR1 and AV14BR2 models, respectively. The AV14UR and AV18UR model predictions are 
almost coincident. The squares (triangles) are the experimental values of ref. || (ref. |10|). The 
straight lines are linear fits of the theoretical results. 



11 



TABLES 
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TABLE I. Singlet a s and triplet at S-wave scattering lengths (fm) for n — 3 H zero energy 
scattering calculated with the AV14 potential and the Kohn (rows labelled K) or Hulthen (row 
labelled H) variational methods. N c is the number of channels included in the CHH expansion 
of the wave functions (the case N c = corresponds to including in the WF only the asymptotic 



terms). The last row reports the results obtained in ref. |21] by solving the FY equations. 
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TABLE II. Total cross section ax (b) and coherent scattering length (fm) for n — 3 H zero energy 
scattering calculated with the AV14UR and AV18UR potential models. The last rows report the 
experimental values. 
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TABLE III. 3 He binding energy B% (MeV) and singlet a s and triplet at s-wave scattering 
lengths (fm) for p — 3 He scattering calculated with the AV18 and AV18UR potential models. The 
last rows report the experimental values. 
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